skip to main content


Search for: All records

Creators/Authors contains: "Heelis, Roderick"

Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher. Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?

Some links on this page may take you to non-federal websites. Their policies may differ from this site.

  1. Daedalus MASE (Mission Assessment through Simulation Exercise) is an open-source package of scientific analysis tools aimed at research in the Lower Thermosphere-Ionosphere (LTI). It was created with the purpose to assess the performance and demonstrate closure of the mission objectives of Daedalus, a mission concept targeting to performin-situmeasurements in the LTI. However, through its successful usage as a mission-simulator toolset, Daedalus MASE has evolved to encompass numerous capabilities related to LTI science and modeling. Inputs are geophysical observables in the LTI, which can be obtained either throughin-situmeasurements from spacecraft and rockets, or through Global Circulation Models (GCM). These include ion, neutral and electron densities, ion and neutral composition, ion, electron and neutral temperatures, ion drifts, neutral winds, electric field, and magnetic field. In the examples presented, these geophysical observables are obtained through NCAR’s Thermosphere-Ionosphere-Electrodynamics General Circulation Model. Capabilities of Daedalus MASE include: 1) Calculations of products that are derived from the above geophysical observables, such as Joule heating, energy transfer rates between species, electrical currents, electrical conductivity, ion-neutral collision frequencies between all combinations of species, as well as height-integrations of derived products. 2) Calculation and cross-comparison of collision frequencies and estimates of the effect of using different models of collision frequencies into derived products. 3) Calculation of the uncertainties of derived products based on the uncertainties of the geophysical observables, due to instrument errors or to uncertainties in measurement techniques. 4) Routines for the along-orbit interpolation within gridded datasets of GCMs. 5) Routines for the calculation of the global coverage of anin situmission in regions of interest and for various conditions of solar and geomagnetic activity. 6) Calculations of the statistical significance of obtaining the primary and derived products throughout anin situmission’s lifetime. 7) Routines for the visualization of 3D datasets of GCMs and of measurements along orbit. Daedalus MASE code is accompanied by a set of Jupyter Notebooks, incorporating all required theory, references, codes and plotting in a user-friendly environment. Daedalus MASE is developed and maintained at the Department for Electrical and Computer Engineering of the Democritus University of Thrace, with key contributions from several partner institutions.

     
    more » « less
  2. Abstract

    Techniques developed in the past few years enable the derivation of high‐resolution regional ion convection and particle precipitation patterns from the Super Dual Auroral Radar Network (SuperDARN) and Time History of Events and Macroscale Interactions during Substorms All‐Sky Imager (ASI) observations, respectively. For the first time in this study, a global ionosphere‐thermosphere model (GITM) is driven by such high‐resolution patterns to simulate the I‐T response to the multi‐scale geomagnetic forcing during a real event. Specifically, GITM simulations have been conducted for the 26 March 2014 event with different ways to specify the high‐latitude forcing, including empirical models, high‐resolution SuperDARN convection patterns, and high‐resolution ASI particle precipitation maps. Multi‐scale ion convection forcing estimated from high‐resolution SuperDARN observations is found to have a very strong meso‐scale component. Multi‐scale convection forcing increases the regional Joule heating (integrated over the high‐resolution SuperDARN observation domain) by ∼30% on average, which is mostly contributed by the meso‐scale component. Meso‐scale electron precipitation derived from ASI measurements contributes on average about 30% to the total electron energy flux, and its impact on the I‐T system is comparable to the meso‐scale convection forcing estimated from SuperDARN observations. Both meso‐scale convection and precipitation forcing are found to enhance ionospheric and thermospheric disturbances with prominent structures and magnitudes of a few tens of meters per second in the horizontal neutral winds at 270 km and a few percent in the neutral density at 400 km through comparisons between simulations driven by the original and smoothed high‐resolution forcing patterns.

     
    more » « less
  3. Abstract

    We report on a new method to derive the on‐orbit electron density using the Tri Global Navigation Satellite System (GNSS) Radio‐occultation System (Tri‐GNSS Radio occultation System (TGRS)) differential total electron content data and compare it to the Constellation Observing System for Meteorology Ionosphere and Climate‐2 Ion Velocity Meter (IVM) ion density data. We found that the IVM ion density is about 8%–15% lower than the TGRS derived density at the insertion orbit (∼710 km) and 5% higher at the mission operation orbit (∼540 km) for reasons that are currently unknown. Using a linear coefficient, we scaled the IVM data to remove the offset between TGRS‐derived electron density and the IVM ion density for the two orbital heights. We believe the scaled IVM densities eliminate any inter‐spacecraft discrepancy, making the IVM data suitable for use in high precision multi‐satellite scientific investigations of longitudinal and local time variations of non‐migrating tides, planetary waves and space weather operational applications.

     
    more » « less
  4. Abstract

    Meso‐scale plasma convection and particle precipitation could be significant momentum and energy sources for the ionosphere‐thermosphere (I‐T) system. Following our previous work on the I‐T response to a typical midnight flow burst, flow bursts with different characteristics (lifetime, size, and speed) have been examined systematically with Global Ionosphere‐Thermosphere Model (GITM) simulations in this study. Differences between simulations with and without additional flow bursts are used to illustrate the impact of flow bursts on the I‐T system. The neutral density perturbation due to a flow burst increases with the lifetime, size, and flow speed of the flow burst. It was found that the neutral density perturbation is most sensitive to the size of a flow burst, increasing from ∼0.3% to ∼1.3% when the size changes from 80 to 200 km. A westward‐eastward asymmetry has been identified in neutral density, wind, and temperature perturbations, which may be due to the changing of the forcing location in geographic coordinates and the asymmetrical background state of the I‐T system. In addition to midnight flow bursts, simulations with flow bursts centered at noon, dawn, and dusk have also been carried out. A flow burst centered at noon (12.0 Local Time [LT], 73°N) produces the weakest perturbation, and a flow burst centered at dusk (18.0 LT, 71°N) produces the strongest. Single‐cell and two‐cell flow bursts induce very similar neutral density perturbation patterns.

     
    more » « less
  5. Abstract. The lower-thermosphere–ionosphere (LTI) system consists of the upper atmosphere and the lower part of the ionosphere and as such comprises a complex system coupled to both the atmosphere below and space above. The atmospheric part of the LTI is dominated by laws of continuum fluid dynamics and chemistry, while the ionosphere is a plasma system controlled by electromagnetic forces driven by the magnetosphere, the solar wind, as well as the wind dynamo. The LTI is hence a domain controlled by many different physical processes. However, systematic in situ measurements within this region are severely lacking, although the LTI is located only 80 to 200 km above the surface of our planet. This paper reviews the current state of the art in measuring the LTI, either in situ or by several different remote-sensing methods. We begin by outlining the open questions within the LTI requiring high-quality in situ measurements, before reviewing directly observable parameters and their most important derivatives. The motivation for this review has arisen from the recent retention of the Daedalus mission as one among three competing mission candidates within the European Space Agency (ESA) Earth Explorer 10 Programme. However, this paper intends to cover the LTI parameters such that it can be used as a background scientific reference for any mission targeting in situ observations of the LTI. 
    more » « less
  6. Abstract

    The Communication/Navigation Outage Forecast System (C/NOFS) satellite's VEFI payload performed frequent recordings of the vector electric field in the band 0–16 kHz during the epoch 2008–2014. The Vector Electric Field Instrument (VEFI) was supported by ion‐composition data from the Coupled Ion Neutral Dynamics Investigation (CINDI) instrument. We focus here on statistics of these “burst‐mode” recordings, of which 6,890 (mostly ~12‐s duration) records meet stringent quality‐control criteria, allowing inference of the wave vectorkand its orientation relative to the Earth's magnetic fieldB0. The 6,890 records occur between ±13° (geographic) latitude and between ~ 400‐ and 850‐km altitude, mostly in the topside ionosphere. The wave activity is dominated by terrestrial lightning. We analyze the whistler‐wave intensity and polarization for each pixel in the time‐frequency spectrogram for each record. We then gather weighted statistics on wave polarization, naturally weighted by wave intensity. In this manner we arrive at statistical results that represent the bulk of the energy flow due to whistler waves. Despite rather nonstationary statistics, we can reach three empirical results.

    We see no evidence of a low‐latitude suppression of whistler‐wave activity, in contrast to the predictions of models of transmission through a laminar ionosphere.

    The wave vector polar angle is always in the range 40° to 90° from parallel toB0. This indicates that the propagation at low latitudes is dominated by oblique, not ducted, whistlers.

    At the lowest magnetic latitudes, the wave vector polar angle with respect toB0becomes nearly 90°.

     
    more » « less